A large Nc perspective on the QCD phase diagram 



D. Toublani 

''Physics Department, University of Illinois at Urbana- Champaign, Urbana, IL 61801 

(Dated: February 1, 2008) 

The transition between the hadronic phase and the quark gluon plasma phase at nonzero tem- 
perature and quark chemical potentials is studied within the large- A^c expansion of QCD. 



I. INTRODUCTION 
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The study of QCD at nonzero temperature and quark chemical potentials is of crucial importance to understand 
a wide range of different physical phenomena, from heavy-ion-collision experiments to neutron stars and cosmology. 
This has led to numerous theoretical investigations of the phase diagram of QCD at nonzero temperature and quark- 
' chemical potentials. 

Historically, numerical simulations using Monte Carlo techniques have been very fruitful for the study of QCD at 
nonzero temperature and zero quark chemical potentials. However, this approach faces a major problem at non-zero 
quark chemical potentials: In general, the fermion determinant is complex and Monte Carlo techniques do not work. 
' This is the so-called 'sign problem'. There is one exceptional case where there is no sign problem: QCD at non-zero 

i—^ isospin chemical potential and zero baryon and strangeness chemical potentials. In this case standard Monte Carlo 
I . techniques can been used. Lattice simulations show that the QCD phase diagram in this particular case is very rich 
Q Q' sign problem represents a severe challenge for our general understanding of the QCD phase diagram 

^SJ . at nonzero temperature and quark chemical potentials. Eventhough there is no general method to solve the sign 
J> • problem, recent advances have been made to circumvent it at small quark chemical potentials [1,^ IE 111' These 
On \ recent studies have concentrated on the transition between the hadronic phase and the quark-gluon-plasma phase 
at nonzero baryon chemical potential and zero isospin and strangeness chemical potentials, and in particular on the 
corresponding critical temperature as a function of quark chemical potentials, /x^, /^s)- 

^~») . At small chemical potentials, two remarkable properties of the critical temperature seem to emerge from lattice 
' studies. The critical temperature weakly depends on the chemical potentials, and Tc{ii,fi,0) ~ Tc(^,— ^,0) Q,|3l- 
These properties are rather puzzling, the second one in particular, since physics at nonzero baryon chemical potential 
' and zero isospin chemical potential is rather different from physics at zero baryon chemical potential and nonzero 
I isospin chemical potential. In this article we shall use the large- iVc expansion of QCD to shed light on these properties 
ft) • and show that they naturally emerge in this context. We shall also use the large- A^c expansion to show that if there is 
[ a first order phase transition between the hadronic phase and the quark-gluon-plasma phase at a critical temperature 
, Tc{fi, —fi,fis) ^ 140 MeV, then there is also a first order phase transition at a critical temperature Tc(/z, /z, /Zs), and 
' vice versa. This is important since simulations at /i„ = —fJ-d can be readily performed using standard Monte Carlo 
techniques, whereas simulations at /i„ = fid suffer from the sign problem. 
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II. LARGE-iV, EXPANSION 



Since its inception, the large- A^c expansion in QCD has been a useful tool both conceptually and phenomenologically 
1^ . It is insightful because it provides us with a small parameter expansion that distinguishes classes of Feynman 
diagrams. Although the theory has not yet been solved at leading order in 1/Nc in four dimension, properties of these 
classes of Feynman diagrams can be used to determine relations between some observables. In this article, we shall use 
the usual large- A"c expansion |^ augmented by recent developments on the large- A'c expansion at nonzero temperature 
and chemical potential In the large- TVc expansion, the leading-order contribution to the pressure is given by planar 
diagrams that contain only giuons, and the next-to-lcading order contribution is given by planar diagrams with only 
one quark loop as a boundary. This expansion is still valid at nonzero temperature and quark chemical potentials, 
provided that the temperature and the quark chemical potentials are of order A'^, and if the different superfluid and 
superconducting phases are avoided 0, 0| . Following the arguments developed for two flavors in |§| , the large- A^c 
expansion leads to the following expression for the pressure in QCD with Nf quark flavors, /, in a finite volume V: 

p{T, {ruf, fif}) = eo({m/, /!/}) + y\n Z{T, {ruf, /i/}) (1) 
p(T, {m/, /i/}) = {po{T) + E/ PI {T, mf,fi}) + 0{^)) , (2) 
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where Z{T, {m/, /Lt/}) is the QCD grand canonical partition function at nonzero temperature, T, for Nf quark flavors 
of mass TO/ and chemical potential /i/, and where eo{{m f , n f}) = —limT^oy^^Z{T,{mf,iJ,f}) is the vacuum energy 
density. In the large- A^c expansion in the equation above, PoiT) contains the leading-order diagrams, which are 
planar diagrams with gluons only, and thus depends on T only. At next-to-leading order there is no mixing between 
the different quark flavors. The diagrams that contribute to pi(T, to/,^/) are planar diagrams that contain gluons 
and only one quark loop as a boundary, and since the gluons do not mix the quark flavors, the quark chemical 
potentials do not mix at next-to-leading order. Furthermore, since Z(T, {to/, /i/}) = Z{T, {m/, —fJ-f}) because of CP, 
pi (T, nif^ii^j:) is even in /x/ . The mixing between the chemical potentials of the different quark flavors appears only 
at next-to-next-to-leading order. 

The {/(l)-axial anomaly could spoil this property of the diagrammatic expansion of QCD in the large-iVc limit. 
However, as shown in [IL] . the triangle anomaly does not depend on the temperature. It is straightforward to 
generalize the argument developed in pT| to show that the anomaly is also independent of the quark chemical 
potentials. Another way to show that the anomaly does not mix the different quark chemical potentials is to use 
chiral perturbation theory in the hadronic pha se at nonzero temperature, zero baryon chemical potential, nonzero 
isospin and strangeness chemical potentials \13i ll4L llSl llq . Chiral perturbation theory is a low-energy effective 
theory of QCD solely based on the symmetry of QCD. It contains only the Goldstone modes due to the spontaneous 
breaking of chiral symmetry: tt, ry, and K. In the large- A'^c limit, the chiral symmetry breaking is different, and the 
rj' becomes a Goldstone boson as well. The large- iVc expansion and the inclusion of the ry' have been implemented in 
chiral perturbation theory .12., JJj . The pressure of the hadronic phase can easily be calculated following ^] . In the 
hadronic phase, we find that the temperature and quark chemical potentials enter the pressure and mix at 0{N^) in 
chiral perturbation theory. Therefore, in chiral perturbation theory the quark chemical potentials mix only at order 
7V° in the large- TVc expansion in the hadronic phase, in complete agreement with the diagrammatic argument sketched 
above and developed in detail in ,9]. We thus conclude that the anomaly does not spoil the diagrammatic argument 
presented in 0. 



III. CRITICAL TEMPERATURE 



At finite volume V, the separation between the hadronic phase and the quark gluon plasma phase corresponds to a 
peak in the specific heat Cy — de/dT\v, where e = —p + Tdp/dT+J2f f^fdp/dfif is the energy density. In the infinite 
volume limit the peak value of Cy stays finite if there is a crossover, and it diverges if the phase transition is first 
order or second order. The scaling of the peak value with the volume depends on the nature of the phase transition. 
Therefore, the critical temperature as a function of the quark chemical potential, Tc(/i„, /x^, /is), is implicitly given by 



dT 



0. (3) 



From (0) and ^ and at zero chemical potentials, we find that the difference between the critical temperature for very 
massive quarks, i.e. pure Yang-Mills theory, Tj^, and the critical temperature for QCD, T^'^^, should be of order 
1/A^c: 

rrQCD _ rpYM i 

' Tj^' =^(Ar)- (4) 

This is indeed what has been observed on the lattice for Nc = 3: Tc — 270 MeV for pure Yang-Mills 0|, and 
Tc ~ 175 MeV for QCD with three flavors [T^ ^2X\l l2ll |. Similarly, the critical temperatures for pure Yang-Mills 
theories with different number of colors have been computed on the lattice, and they have been found to differ by 
0{\/N'^), in agreement with the large- iVc expansion, 

T^^{N,) , 0.76(6) ^ 



rpYM ' ]\T2 ' ' 

where is the critical temperature for pure Yang-Mills theory when Nc oo [T8l |. 

At nonzero quark chemical potentials the large-iVc expansion in (jSJ implies that for QCD with iriu = rud, the 
critical temperature as a function of quark chemical potential, Tc{fJ.u, fJ'd, fJ's), must satisfy the following relation 

TciiJ., Us) 
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Therefore, for /Zs = 0, the critical temperature that separates the hadronic phase and the quark-gluon-plasma phase 
at nonzero baryon chemical potential, fj,}, = {jiu + Md)/2, and zero isospin chemical potential, /i^ = — /Xd)/2, differs 
from the critical temperature at zero baryon chemical potential and nonzero isospin chemicalpotential by This 
is in complete agreement with recent results obtained in numerical lattice simulations [ij, |^ , and in various models 
|22, 12,24]. This is a useful relation since lattice simulations at /if, = and /i^ ^ do not suffer from the sign problem 
present at nonzero /ib, and are therefore much easier to perform. 

More explicitly and for QCD with three quark flavors with to,; = ^ nis in a finite volume, we can perform 
a Taylor expansion of dCy/dT around the critical temperature that separates the hadronic phase from the quark- 
gluon-plasma phase at zero chemical potentials, Tg. Using (0), we find that 



dCy 
dT 



= ai 




bo + bi 



T-Tn 



To 



rp2 



bso + b. 



T-Tn 



sl- 



To 



rp2 

^0 



(7) 



Notice that the coefficients related to /is differ from those related to /in and /i^ since m„ = irid mg. Using (O to 
solve (O, we find that the critical temperature as a function of the quark chemical potentials is given by 



1 / bo fil + 



ai 



rp2 

^0 



(8) 



Therefore, we find that the large-A'c expansion leads to interesting insight on the critical temperature as a function 
of quark chemical potentials. First the curvature of the critical temperature for small chemical potential is 1/Nc 
suppressed. Second, for a given number of colors TVc, the curvature of the critical temperature as a function of baryon 
chemical potential should increase with the number of flavors. For Nf degenerate quarks, this increase in the curvature 
should be linear in Nf up to 0{1/Nc) corrections. This has been indeed observed in different lattice simulations for 



2, 3, and 4 



n iVf 

m 



Third, the leading-order dependence in the quark chemical potentials is even in /i„, /i^;, and 



/is separately, and does not mix them. The mixing term appears at next-to-next-to-leading order only. Therefore 
effects due to the mixing of the quark chemical potentials should be of the order of This expression for the 

critical temperature in the large- A'^c expansion agrees with and sirnply explains several results that have been found 
in numerical simulations and in various models Q, Q |^, IE Hji 1^ • 



IV. ORDER OF THE PHASE TRANSITION 



The arguments developed above do not depend on the nature of the phase transition. This is, however, an important 
question to address. Indeed, for QCD with three colors and physical quark masses, the separation between the hadronic 
phase and the quark-gluon plasma phase at nonzero baryon chemical potential, zero isospin and strangeness chemical 
potentials is believed to be a crossover at low chemical potential and a first order phase transition at higher chemical 
potential [2^ [2^ ll^l • Several recent lattice simulations have found the critical endpoint that corresponds to the end 
of this first order phase transition line [^Q. However, these different simulations yield results that significantly differ 
on the precise location of this critical endpoint. 

Lattice simulations at nonzero baryon chemical potential suffer from the sign problem. No such problem is present 
at nonzero isospin chemical potential and zero baryon and strangeness chemical potentials. It is therefore easier to 
perform lattice simulations in the latter case. In the section above, we have shown that the large- A'c expansion leads 
to a relation between the critical temperature at /i„ — fj,d with the critical temperature at /x„ = —fj,d that is valid 
at next-to-leading order. It is therefore natural to investigate if such arguments can lead to a relation between the 
nature of the phase transition in these two cases. 

In the infinite volume limit, if the phase transition is first order, then there is a latent heat, Lh — Tc disc s, where s 
is the entropy density. The latent heat is related to the discontinuity of the quark-antiquark condensate, (gg), through 
the Clausius-Clapeyron relation derived in [2^12^ 



Lh = a,rj. ir.\ disc {qq), (9) 

where Q\x means that quantity Q is evaluated at constant x. Starting from 0, and using the same reasoning that leads 
to ©, we find that dTd d7nq\{^j^ ~ f/-/Vc- At zero temperature the large- TV^expansion leads to {qq)T=o ~ 0{Nc)- 
At nonzero temperature, /x„ — ^fid E^nd /is = 0, based on lattice simulations and on chiral perturbation theory 

|l6l Isof , we know that (qq) is almost independent of the chemical potential and of the temperature in the hadronic 
phase, provided T < 140 MeV. For instance in chiral perturbation theory at zero chemical potentials, it was found 



4 



that for T < 140 MeV, (qq) / {qq)T=o ^ 80% for QCD with two massive quarks, and with corrections due to massive 
non-Goldstone modes taken into account l3C!'|. Thus, from both lattice simulations and chiral perturbation theory, we 
have that (qq) ~ {qq)T=o — 0{Nc) for T < 140 McV in the hadronic phase. Therefore, we conclude that if a first 
order phase transition between the hadronic phase and the quark-gluon plasma phase takes place at T ^ 140 MeV, 
the corresponding latent heat is of order N^. This conclusion cannot be reached for higher temperatures since the 
value of the quark-antiquark condensate decreases and is no longer 0{Nc). 

In the large-A^c perspective described above and if there is a first order phase transition at T < 140 MeV, we 
can use the same reasoning that led to Q to show that the latent heat as a function of quark chemical potentials, 
Lhipu, fJ'd, fJ's), has to satisfy the following relation 

L,i^,,^^,^^.) ~Ui^, ~^,,^^.) ^ 1 forr<140MeV. (10) 

In other words, if there is a first order phase transition between the hadronic phase and the quark-gluon plasma 
phase at a temperature Tc{n, —ii,fis) ^ 140 MeV, with a latent heat = 0{N^), or equivalently disc {qq) = 0{Nc), 
then, according to the large- iV^ expansion, equations 10 and (|10|l . there should be a first order phase transition at 
a temperature Tc(/i, /u, /is) ^ T'c(/i, — /x,/is) -I- 0{1/N^). Thus if lattice simulations at /!„ = —fid were to find a first 
order phase transition at temperatures below ~ 140 MeV, then the l/N^ expansion predicts that there should also 
be a first order phase transition at the same fiu = fJ-d and at the same temperature, up to 0{1/N^) corrections. 



V. CONCLUSIONS 



We have used the large- A'c expansion of QCD to study the phase transition between the hadronic phase and the 
quark-gluon plasma phase at nonzero temperature and quark chemical potentials. We have shown that the critical 
temperature depends on the chemical potentials at next-to-leading order in the large- A^c expansion. We have also 
shown that there are relations between the critical temperature at nonzero baryon chemical potential and zero isospin 
chemical potential, and the critical temperature at nonzero isospin chemical potential and zero baryon chemical 
potential. These relations are valid at next-to-leading order in the large-iVc expansion. Finally, based on large- 
Nc arguments, we have shown that it should be possible to determine relatively accurately the position of a first 
order phase transition between the hadronic phase and the quark-gluon plasma phase in the QCD phase diagram 
at nonzero baryon chemical potential and zero isospin and strangeness chemical potentials, by performing ordinary 
lattice simulations at nonzero isospin chemical potential and zero baryon and strangeness chemical potentials, as long 
as the latent heat is large enough. This should be a valid strategy up to temperatures of ~ 140 MeV where the 
quark-antiquark condensate is not significantly different from its value at zero temperature. 
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